Optimum unambiguous discrimination of two mixed quantum states 
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We investigate generalized measurements, based on positive-operator-valued measures, and von 
Neumann measurements for the unambiguous discrimination of two mixed quantum states that 
occur with given prior probabilities. In particular, we derive the conditions under which the failure 
probability of the measurement can reach its absolute lower bound, proportional to the fidelity of 
the states. The optimum measurement strategy yielding the fidelity bound of the failure probability 
is explicitly determined for a number of cases. One example involves two density operators of rank 
d that jointly span a 2d-dimensional Hilbert space and are related in a special way. We also present 
an application of the results to the problem of unambiguous quantum state comparison, generalizing 
the optimum strategy for arbitrary prior probabilities of the states. 
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Many applications in quantum communication and 
quantum cryptography are based on transmitting quan- 
tum systems that, with given prior probabilities, are pre- 
pared in one from a set of known mutually nonorthogonal 
states. Since perfect discrimination between nonorthog- 
onal quantum states is impossible, measurement strate- 
gies for state discrimination have been developed that are 
optimized with respect to various criteria Q. Here we 
consider unambiguous discrimination, requiring that the 
outcome of the measurement be error-free. For two mixed 
quantum states unambiguous discrimination is possible 
with a finite probability of success if the supports Q 
of their density operators are not identical. When the 
measurement fails, it returns an inconclusive answer but 
never an error. In the optimal measurement strategy the 
failure probability is minimum. 

The problem of unambiguously discriminating mixed 
quantum states arises for instance when given pure states 
undergo a specified decoherence process during transmis- 
sion through a quantum channel, or when the quantum 
system is known to be in a pure state which has to be 
assigned to a particular set out of a number of given 
sets of pure states, each set corresponding to a mixed 
state. While for two pure states the minimum failure 
probability has been known since long 0, Q , the study 
of unambiguous discrimination among mixed states, or 
sets of pure states, respectively, started only recently 
HH 0,110 El Ell A complete solution, determining 
the minimum achievable failure probability for arbitrary 
prior probabilities of the states, has been obtained for the 
special cases of discriminating a pure and a mixed state 
[a, and of two mixed states of rank d in a (d + 1)- 
dimensional joint Hilbert space 0. For discriminating 
two arbitrary mixed states, bounds have been derived for 
the failure probability yj , in terms of the fidelity of the 
states. In this paper we perform a more detailed analy- 
sis, investigating the conditions under which the lowest 
bound, proportional to the fidelity, can be reached, and 
deriving also the von Neumann measurements for unam- 
biguous discrimination. 



We start by recalling that a measurement for distin- 
guishing two quantum states, characterized by the den- 
sity operators pi and pi and the prior probabilities r\\ 
and 772 = 1 — ^71 , respectively, can be formally described 
by three positive operators with J2k=o ^fe = I> where 
/ is the identity. These detection operators are defined in 
such a way that Tr^pH^) with fc = 1,2 is the probability 
that a system prepared in a state p is inferred to be in 
the state pk, while Tr(pHo) is the probability that the 
measurement fails to give a definite answer. When all 
detection operators are projectors, the measurement is a 
von Neumann measurement, otherwise it is a generalized 
measurement based on a positive operator-valued mea- 
sure (POVM). From the detection operators schemes 
for realizing the measurement can be obtained |l2t Il3| . 

It is our aim to investigate the optimum measurement 
strategy that minimizes the total failure probability 

Q = 77iTr( Pl no) + 77 2 Tr(p 2 n ). (1) 

From the relation between the arithmetic and the ge- 
ometric mean and from the Cauchy-Schwarz-incquality 
[HQ it follows that Q > 2 v / T ll T 1 2Tiip 1 U a )TT{p 2 U a ) > 
2 v / 77r^2 Maxy \Tv{U ^/p\\Iq^/p2)\, where U describes an 
arbitrary unitary transformation. The failure proba- 
bility takes its absolute minimum when the two equal- 
ity signs hold. This is true if and only if both the 
relations 77iTr(p 1 n ) = ^Tr^Ho) and U^fplyJTl^ ~ 
y/piVHo are fulfilled. From the first relation we con- 
clude that the number of inconclusive results is equally 
distributed among the two incoming states. After mul- 
tiplying the second relation with its Hermitean conju- 
gate, the two conditions for equality can be combined to 
yield v / Ho(»72P2 - ViPi)V^o = 0. Since in the POVM- 
formalism the detection operators transform a quantum 
state according to p — + ^ fe a/TTaT p VHfc ,13], it follows 
that the total failure probability is smallest when in case 
of failure the two density operators are transformed into 
states that are identical after normalization and therefore 
cannot be further discriminated. 



2 



We now recall that unambiguous discrimination of two 
states leads to the requirement piH 2 = p 2 Hi = 00. 
Substituting n = / — III — ±1 2 into the inequality for the 
failure probability Q |lCj . given above, we arrive at 

Q > 2^nrnM&xu \Tr(Uy/pI^)\ = F, (2) 

where F = Tr [(y^ pi^) 1/2 } is the fidelity 0. Using 
a different method, it has been found already previously 
by Rudolph et al. that 



{ Vmin + VmaxF 2 otherwise, 



(3) 



with ?/ m i n (?7 m ax) denoting the smaller (larger) of the prior 
probabilities. Here, in addition, we obtained the neces- 
sary and sufficient conditions that the detection opera- 
tors have to fulfill in order to reach the fidelity bound 
Qo- They can be summarized as 

n = / - n : - n 2 > o, ri! > o, n 2 > o, (4) 
Pl u 2 = P2 n 1 = o, (5) 

r?iTr( Pl n ) = rh[l - Tr(piIIi)] = y^F, (6) 
77 2 Tr(p 2 n ) = r? 2 [l - Tr(p 2 II 2 )] = ^/rjrq^F. (7) 

In the following we investigate the conditions under 
which detection operators exist that satisfy Eqs. Q - 
(0) . For this purpose we use the spectral representations 



di 

Pi = E r; l r ^ n l' 
i=i 



P2 



d 2 

E 

m=l 



where n,s m ^ 0, and (n|r m ) = 5i, m = (si\s„. 
more, we introduce the projection operators 



(8) 
Further- 



P^ = 



di 

E 

i=i 



n)(n\, p 2 = Y, \sm)(s m \, (9) 

m= 1 



where pi\fj) = 0. The states {(r^)} form an orthonor- 
mal basis in the ^-dimensional subspace of 7iio that is 
spanned by all states that are orthogonal to Pi, where 
d\ = di\i + di± — d\. The identity is then given by 



I = Pi. 



P?=P 



All 



A 



P 1 +P 1 +P! 1 . (12) 



Here the operator P 2 = I — Pu_ — Pn projects onto the 
subspace TL' 2 spanned by those states that are orthogonal 
to both Pu_ and Piy, implying that piP 2 = 0. Instead of 
decomposing the eigenstates of p\ , we might as well have 



started from Is, 



Pik 



|s^j), obtaining instead of 



Eq. (|12JI the alternative decomposition 
/ = P 2 ± + P = P 2 ± + P 2 || + P[ = P 2 + P 2 + P{, (13) 

where the projectors are defined analogously. 

Now we can specify the general structure of all detec- 
tion operators, Hi and n 2 , that describe unambiguous 
discrimination, i. e. satisfy Eqs. (0J and JSJ. We write 



ni = E«»^l 



where < a' < 1 and 



di± 

E a *. 



\Vi){Vj 



(14) 
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with {uji} being 



a unitary matrix. We note that J^j \ v j)( v j\ = A_L smce 
the eigenstates \v'j) form a complete orthonormal basis 
in 7ii_\_. For representing n 2 we start from the same 
decomposition of the identity, and take into account that 
none of the eigenstates of Ho must be contained in the 
subspace H' 2 when the failure probability is to be as small 
as possible. This leads to 



di di 

n 2 = Ymm+pi = E /wfo 



-/-P11 



(15) 



i=i 



and the nonnormalized states |r! ) = P 2 |n)- We can con- 
struct a complete orthonormal basis {\hk}} in the sub- 
space 7ii\\ spanned by the state vectors P 2 |n), using the 
recursion relation \h k ) = P 2 \r k ) - J2i=i \hi){hi\P 2 \r k ) 
and determining \h k ) — \h k ) /\\hk\\ Q- The dimen- 
sionality d\\\ of H\\\ is equal to the rank of the matrix 
formed by the elements (ri\P 2 \r n ). Similarly, in the sub- 
space 7iij_ that is spanned by the nonnormalized vectors 
\r\ ] ) = (/ — P 2 )|n), we can obtain an orthonormal ba- 
sis of dimension di±. The respective projection 
operators into the two orthogonal subspaces are 

din d 1± 

Pi\\=J2\h k ){hk\, Pi±=^2\vi)(vi\, (10) 
fc=i t=i 

where p 2 \vi) = 0. The operator P w = Pn + Pi± projects 
onto a subspace TiiQ of dimension di\\ + d%±. Noticing 
that Tr[(Pi — Pi)pi] = 0, we construct the operator 



Pi = Pn 



P 



1_L 



Pi 



E 



r j)( r j\ 



(11) 



where < # < 1 and J2i=i = A- The constants 

ctij and 0ij are subject to the constraint that Ho > 0. 
Clearly, when 

Pi = -Pill = I, and consequently also P 2 = P 2 ii = I, 
it follows that Hi = n 2 = and Ho = /, yielding a 
unit failure probability that makes error-free discrimina- 
tion impossible. We therefore require that Pi_l 0, or 
P 2 ± 7^ 0, respectively, which, because of normalization, 
is equivalent to 

Tr(P lP2 ) < Tv(P n p 2 ), Tr(P 2Pl ) < Tr(P 2 || Pl ). (16) 

Before studying the optimum measurement, let us con- 
sider the von Neumann measurements for unambiguous 
discrimination. If a'j = for all j, and P[ = 1 for all i, it 
follows that ni = and n 2 = A + P£. Hence n = P u 
with the failure probability Qni = + ?/ 2 Tr(Pip 2 ). 
Another von Neumann measurement is generated when 
a'j = 1 for all j, and j3[ = for all i, giving Hi = Pi± and 
n 2 = P 2 . Then Hq = Pin, with the failure probability 



Q 



Will 



7?iTr(P 2 pi) + ?7 2 Tr(P 1 ||p 2 ), 



(17) 



3 



where the relation Tr(Pi|| (Oi ) = l-Tr(P 1± pi) = Tr(P 2 pi) 
has been applied. In this measurement the state is un- 
ambiguously found to be pi when a detector click occurs 
in a direction orthogonal to all eigenstates of pi. On the 
other hand, for a click in a direction orthogonal to both 
Piy and Pu_, the state is determined to be p2 with cer- 
tainty, and in the rest of cases the result is inconclusive. 
So far we relied on Eq. 1|12|) . Based on the complemen- 
tary decomposition of the identity, Eq. (|13(l , we obtain an 
alternative pair of von Neumann measurements. These 
yield the failure probabilities Qn2 = V2 + 7 7iTr(P 2 pi) and 



QjV2|| = »72Tr(Pi/92) + 77iTr(P 2 ||/9i) 



(18) 



Obviously Qjv2|| < Qni and Qjvi|| < Qn2- 

We now return to the optimum measurement. Since 
the von Neumann measurements can be performed for 
arbitrary given parameters, the optimized failure proba- 
bility certainly obeys the inequality 



Qopt < Min{Qivi||,QiV2||}. 



(19) 



According to Eqs. and the absolute minimum of 
the failure probability, Qo = 2^/r\\r\i F, is reached if and 
only if the two conditions Tr(piIIo)/P = ^772/771 and 
P/Tr(p 2 n ) = \Jr\ijrii are fulfilled. However, due to the 
structure of the operators 111 and II2 , the possible values 
of Tr(p fc n ) = 1 - Tr(p fc n fc ), for jfe = 1,2, have a lower 
bound. In particular, 

Tr( Pl n ) > 1 - Tr(P 1±Pl ) = Tr(P 2 pi), (20) 
Tr(p 2 n ) > Tr(Pi||p 2 )-Tr(Pip 2 ) =Tr(Pip 2 ),(21) 

where in the first equation the equality sign holds when 
a'j = 1 in Eq. I|14l) . and in the second equation the 
equality is reached when (3[ = 1 in Eq. JTSJ. Therefore 
we obtain that the condition, 



Tr(P 2( 0i) 



< 



< 



Tr(P lP2 ) 



(22) 



is necessary, i. e. the fidelity bound, Q = Qo, can only 
be reached in part or in the whole of this interval. 

The interval specified by Eq. ll'L'l) is not empty only 
when Tr(P 2/ 9i)Tr(Pip 2 ) < F 2 . For two density opera- 
tors that violate this inequality, the failure probability 
Qo cannot be achieved for any values of the prior prob- 
abilities of the states, and the conditions and are 
then of no help for determining the optimum measure- 
ment. Moreover, our result shows that in general the 
lower bound Qo can only be reached in an interval of the 
ratio 772/771 that is smaller than the interval given in Eq. 
©, since Tr(P 2 pi)/P > F and P/Tr(Pip 2 ) < 1/F. The 
latter relations follow from the general inequalities 

Tr(P 2 pi)Tr(Pi|| P2 ) > F 2 , Tr(Pip 2 )Tr(P 2 ||pi) > F 2 

(23) 

that can be readily inferred from Eqs. (|17fl . I|18|) and 0. 

The parameter intervals in Eqs. and (1221 co- 
incide when Tr(Pip 2 ) = Tr(P 2 pi) = F 2 . This condi- 
tion is fulfilled, e. g. for density operators of the 



form pi = n\n)(ri\ and p 2 = £\ =1 n\s l )(s l \, with 

(ri\sj) — bdij, where the corresponding eigenvalues are 
identical. The fidelity is then found to be F = \b\. 

Another simplification arises when Phi and P 2 )| are 
one-dimensional projectors, d\\ \ = d 2 y = 1. In this 
case equality holds in Eqs. (|23() fla ] which implies that 
F 2 = Tr(P 2 pi)Tr(Pi||p 2 ) > Tr(P 2 pi)Tr(Pip 2 ), where Eq. 
(|16|l has been taken into account. Hence again for any 
two density operators the necessary condition l|2^|l is ful- 
filled for a certain range of the ratio ?7 2 /771 . At the lower 

limit of this range, i. e. for . j^- = Tr ^ lP2 - > , we can write 

Z^JrfJteF = Vi Tr^p^ + %Tr(Pip 2 ) = Q N1 \\, and simi- 
larly we find that at the upper limit 2^/771 772 F = Qn2\\ ■ 
Thus, if Q = Qo in the entire range in Eq. (|22|l . the com- 
plete solution for the optimum measurement is known. 

In general, in order to find the optimum measurement 
strategy that yields the failure probability Qq, we have 
to determine the parameters onj and /3y in Eqs. I|14(l and 
(|15H that satisfy the necessary and sufficient conditions 
0-0. In the following we apply this method to a 
number of special cases. 

First we consider two density operators of rank d in 
a 2c?-dimensional joint Hilbert space. In such a case 
P 2 = and the identity can be alternatively expressed as 
/ = P1+P1 or I = Pu_ + P 2 which means that Pi|| = P 2 , 
P 2 || = Pi and Pi = P 2 _l. We start from Eqs. ® with 
d\ = <i 2 = d and assume that \si) = (|rj) + \fi))/^/2, and 
\vi) = (In) - \n))/V2 {i = l,...,d). Then we obtain 
F = E» and Tr(Pip 2 ) = Tr(P 2 pi) = 1/2. It is 

important to note that in general there exist sets of eigen- 
values {r-j} and {si} where F 2 < 1/4 and the necessary 
condition, Eq. Il22[| . cannot be fulfilled. In the following, 
however, we restrict ourselves to the special case that 
Ti = Si for i = 1, . . . , d, for which F = 1/ \/2. The neces- 
sary condition for the lower bound Qo to be achievable 



then reads 



s/2 



< 



tions ctij — a Sij 



< Further, we find the solu- 



and Pij — (3 5^ (i,j = 1, . . . , d), where 



= and/3 = 2-. 



no has two eigenvalues, 
A = and A x = 2 — a — (3, each with a d fold degener- 
acy. Thus the optimum n is always an operator of rank 
d. Note that 2^/2 — 2 < Ai < 1 in the whole interval 
F < w— < i. Hence in this parameter interval the op- 



t) 2 — f- 

timum detection operators yielding the lower bound Qo 
are Hi = aPu_ and n 2 = /3P 2 j_ ■ At the upper and lower 
limits of the interval the measurement turns into the von 
Neumann measurements that give the failure probabili- 
ties Qjvill = Qn2 and Qat 2 || — Qni, respectively. Since 
in our example Tr(P 2 pi) = Tr(Pip 2 ) = F 2 , we find that 
Qni = 1I1+V2F 2 and Qn2 = V2 + ViF 2 - Thus we derived 
a measurement strategy that yields the equality sign in 
Eq. for two mixed states. 

In our next examples we focus on the case din = 



'■2 1| 



1. First we assume that the density operators 



given in Eq. have arbitrary ranks d\ and d 2 , and 
that (ri\s m ) = aSi^Sm^i with \a\ < 1. This yields 



4 



F = y/^n\a\, Tr(P 2j0l ) = F 2 / Sll Tr(P lP2 ) = F 2 / n 
and dij\— g?2|| =1- For the parameter range specified in 
Eq. (1221) we obtain the optimum detection operators 



IIi 



m F \ \vi)(vi\ 



n 2 



rym; (1-M 2 ) 2 
IWil 



;=2 



d 2 



772 V (1 - |a| 2 ) 2 

' x ii/ rn—2 



El r ')^l ( 24 ) 
(25) 



3 771/ \°?7l | J 



Fi 



•olsi) and IfO = Isi) - 



where we introduced 

This solution can be applied to the problem of quan- 
tum state comparison where two identical quan- 
tum objects are each prepared either in the state \ipi), 
or in the state \ip2), and where we wish to determine 
unambiguously whether the states are equal or differ- 
ent. The task amounts to distinguishing the two-particle 
states pi = \ ViXV'i^il + 1^2, ^2)^2, V"2 1) and 
P2 = ^Qipi,ip2){ipi,ip2\ + \ip2,i>i){il>2,ipi\), where F = 
|(V'i|V'2)|- Upon determining the eigenstates, we find that 
the structure of p\ and pi corresponds to the one treated 
in the above special example, with T\ = s\ = (1 + F 2 )/2 
and \a\ = 2F/(1 + F 2 ). The minimum failure probability 
in unambiguous quantum state comparison follows to be 



Qo 



pt 



2^m^F 

2F 2 



+ Vn 



if 

otherwise 



> 



2F 

T+W 



(26) 

Here rj m i n (77 max ) is the smaller (larger) of the values 771 = 
p\ + p\ and 772 = 2pip2, where p\ and p^ are the prior 
probabilities of the states l^i) and \tp2), respectively. 



As our final example we mention the problem of dis- 
criminating a pure state, p\ = from a mixed 
state P2, or from a set of pure states, respectively, that 
has been introduced as quantum state filtering [5J, Il7j . In 

this case Tr(Pip 2 ) = F 2 and Tr(P 2 pi) = ||r||| 2 . In the 
parameter interval given by Eq. I|22(l the optimum detec- 
tion operators take the form IL = (1 — , f^-F) l" 1 ^" 1 ! 

V V vi J i-ik!n 2 



and n 2 = 1 - 



1-lkfll 2 



P 2 , and the previ- 



ous solution for the minimum failure probability in opti- 
mum unambiguous quantum state filtering 0, [fij is read- 
ily regained. 

In summary, we performed a detailed analysis of the 
probabilistic measurement for unambiguous discrimina- 
tion between two arbitrary mixed quantum states. We 
derived general analytical relations that depend on five 
quantities characterizing the mutual relationship of the 
density operators of the states. These quantities are the 
expressions Tr(Pip 2 ) and Tr(P 1 ||p 2 ), as well as Tr(P 2 pi) 
and Tr(P 2 ||/Oi) and, most importantly, the fidelity F. We 
also showed that the method developed in this paper can 
be used to find complete analytical solutions that de- 
scribe the optimum measurement for special cases. 
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